Self-Focussing Dynamics in Monopolarly Charged Suspensions 



O 

o 

(N 

> 
O 

o 



I 

O 

o 



> 

o 
o 

-I— > 



I 

o 
o 



X 
S3 



Stephan M. Dammer and Dietrich E. Wolf 
Institut fiir Physik, Universitdt Duisburg-Essen, 4'^048 Dmsburg, Germany 
(Dated: February 2, 2008) 

The Smoluchowski equation for irreversible aggregation in suspensions of equally charged particles 
is studied. Accumulation of charges during the aggregation process leads to a crossover from power 
law to sub-logarithmic cluster growth at a characteristic time and cluster size. For larger times the 
suspension is usually called stable, although aggregation still proceeds slowly. In this regime the 
size distribution evolves towards a universal scaling form, independent of any parameter included 
in the theory. The relative width falls off to a universal value o-J?°~0.2017 that is much smaller than 
in the uncharged case. We conjecture that cr^ is a lower bound for the asymptotic relative width 
for all physical systems with irreversible aggregation. 
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In many fields of science irreversible aggregation (or 
coagulation, agglomeration, clustering) phenomena are 
important. Examples include aerosol coalescence, poly- 
merization or gelation (see, e.g., 0, SSSIll and refer- 
ences therein). Such systems have been described using 
Smoluchowski's coagulation equation and extensive an- 
alytical and numerical studies have been carried out for 
uncharged systems 0, [3| • In spite of the fact that charges 
are present in a lot of aggregating systems, much less is 
known for the charged case. In this Letter we investigate 
the influence of charges on the aggregation process. In 
contrast to previous work which studied aggregation be- 
tween oppositely charged particles we investigate the 
aggregation dynamics of equally charged particles with 
the focus on suspensions. A practical example is given in 
where coating of medical drug particles with nano- 
particles is investigated, in order to make them inhalable. 
As an intermediate step it is desirable to keep each pow- 
der component in suspension separately, controlling its 
aggregation by charging the particles equally. They col- 
lide due to Brownian motion, and strong van der Waals 
forces result in irreversible aggregation (Brownian coag- 
ulation). However, the results we derive are valid for 
more general systems, as well. Asymptotically neither 
the solvent nor the diffusion properties are important, as 
long as particles may have high kinetic energy, albeit as 
rare events, and charges accumulate in the aggregation 
process. 

Usually a suspension is regarded as stable with respect 
to the aggregation of equally charged particles, if the par- 
ticle diameter d is small compared to the Bjerrum length 
^B=<Z^/47reoei.fcBT, where q is the charge of the particles, 
T and are the temperature and the relative dielectric 
constant of the carrying fluid, fee is Boltzmann's constant 
and eo is the dielectric constant of vacuum. The reason 
is that for d<£B the energy barrier for bringing two par- 
ticles into contact is larger than the thermal energy. 

Consider the case, where a suspension is initially un- 
stable (as in 0), which means that initially aggregation 
events are frequent. We are going to show that charge 



accumulation eventually leads to a crossover to slow ag- 
gregation at characteristic values of time tc and cluster 
size Sc- Above tc the Bjerrum length becomes larger than 
the aggregate diameter, so that one would regard the sus- 
pension as stable for later times. However, the cluster 
size distribution keeps evolving slowly and approaches 
a universal self-similar form which is independent of any 
parameter included in the theory. A particular important 
consequence of this is that the relative width of the clus- 
ter size distribution starts decreasing at tc and asymp- 
totically reaches a universal value (7^ss0.2017. In this 
sense, the cluster size distribution becomes narrower and 
we refer to this effect as charge-induced self-focussing. 
For Brownian coagulation the value (T^«0.2017 must be 
compared to a^^l for the uncharged case. We conjec- 
ture that (Tj°° is a lower bound for the asymptotic value 
of (Tr for all physical systems with irreversible aggrega- 
tion. Our results are valid as long as aggregation events 
occur frequently compared to dissociation of clusters (ne- 
glected in the following). 

We investigate Smoluchowski's coagulation equation 
for monodisperse initial conditions, i.e. initially all par- 
ticles have the same mass m* and radius a* and carry 
the same charge q* . As mass and charge of aggregates 
are proportional to each other in this case, one index i is 
sufficient in the coagulation equation 
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ni{t) denotes the number density at time t of clusters 
of mass mi=im*, charge qi=iq* and radius ai=i°'a*. 
1/a is the fractal dimension of the aggregates. If they 
are spherical (e.g., coagulating droplets) one has a=l/3. 
More generally one may assume that l>a>l/3 for fluffier 
aggregates. The initial conditions are ni(t=0)=l and 
ni{t=0)=0 for i>l. Rij is the rate coefficient (coagula- 
tion kernel) for mergers between clusters with indices i 
and j. Eq. ^ neglects spatial fluctuations and corre- 
lations. In this sense it is a mean field equation. We 
checked, that a Brownian dynamics simulation of aggre- 
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gating charged spheres in three dimensions is in accor- 
dance with the mean field description |7J. 

The rate coefficients i?^ depend on the particular 
physical system under consideration. In suspensions clus- 
ters of mass rrii diffuse with diffusion constant Di. Usu- 
ally the van der Waals interaction is represented as a 
sticky contact force, leading to irreversible aggregation. 
Therefore, aggregation with a cluster of mass rrij oc- 
curs when the distance between the clusters equals Ui+aj 
(Brownian coagulation). A finite range of the van der 
Waals interaction would basically increase this 'contact 
distance' by a constant, which becomes less and less im- 
portant for growing clusters. For Brownian coagulation 
Rij can be calculated by solving the diffusion equation 
in the presence of an absorbing sphere 0, leading to 
Rij=A'K{ai + aj){Di + Dj)=rij. If Einstein's formula ap- 
plies, one can insert Di~k^T /Qirriai, with the viscosity 
77 of the carrying fluid. This leads to 
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For most uncharged coagulating systems discussed in 
the literature the coagulation kernels arc homogeneous 
functions of degree A, i.e. Rubj—b'^Rij, with A<2 
Eq. for instance, corresponds to A=0. It turns out 
that A is the most important parameter for the uncharged 
case. Depending on its value, the solutions of the co- 
agulation equation are qualitatively different. For A<1 
the cluster size distribution decays exponentially for large 
sizes at all times. In the scaling limit the solutions are 
self-similar and (for A<1) have the form 



n.^{t) = s{ty^(t>{i/s{t)) with s(t) - t~ 



(3) 



The scaling function cj) depends on the details of the un- 
derlying coagulation kernels. The average size s(t) is a 
measure for the average number of primary particles per 
cluster. For 2>A>1 the cluster size distribution develops 
a power-law tail in finite time t at large sizes that vio- 
lates mass conservation {runaway growth)^ which is 
associated with the occurrence of a gelation transition. 

A derivation similar to the uncharged case leads to the 
coagulation kernel for Brownian coagulation of charged 
clusters 
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Wijrij , Wi-j = Kij/ [exp(Kjj) - 1] , 



(4) 



where r.y is the coagulation kernel for uncharged clusters, 
Eq. (O . It was assumed that at contact of the clusters 
the charges are separated by a distance ai-\-aj. 

The important feature of Wij is that collisions of 
strongly charged particles (k^ -^1) are exponentially sup- 
pressed. We are going to show that Wy = exp(— Kij) 
(which is found in charged granular gases |^) leads to 
the same asymptotic results. 



The purpose of this Letter is to study the behav- 
ior of Eq. depending on the parameters A, a and 
k=^J q*'^ / {An€oei-a* k^T) , for rate coefficients that obey 
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exp(Ky) - 1 

or % = (zjy/2 exp(-Ky ) (6) 
k^ij 



with 
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These expressions are motivated by the observations pre- 
sented above. Kij is proportional to the Coulomb energy 
at contact of the clusters divided by a temperature-like 
variable. Eq. |SJ| is the dimensionless coagulation kernel 
for Brownian coagulation of charged clusters, Eq. 
with 2kBT/3r]=l. Whereas this corresponds to A=0 as 
was pointed out above, we investigate the influence of 
Ay^O using the general rate coefficient Eq. ©. 

Let the system initially be unstable (fc^^l) which 
means that in the beginning the exponential terms in 
Eqs. H5lti|l are close to unity. The system basically be- 
haves as if it was uncharged. This changes as soon as 
values of Kij~l become important so that the expo- 
nential terms must be taken into account. Then fur- 
ther aggregation is exponentially suppressed and clus- 
ter growth becomes very slow. This happens when the 
average cluster size s approaches a characteristic value 
Sc given by k^^s^wI. Note that even systems which in 
the uncharged case exhibit runaway growth, i.e. 1<A<2, 

eventually cross over to slow aggregation. According to 

1 

Eq. 10) the crossover time tc is given by Sc^tc^^ for A<1, 
leading to 
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t* for A < 1 



(8) 



where t* is the appropriate time unit. Note that the 
expression for Sc is also valid for A>1. Though similar 
arguments may be applied for tc and A>1 we do not 
discuss this case here. 

The rate equations ^ were solved numerically for var- 
ious choices of the coagulation kernels Rij and different 
values of k. Fig. ^ shows data for Brownian coagulation 
of charged spheres (Eq. ||SJ| with a=l/3). The average 
size s{t) was calculated as s=Mi/Mq with the moments 
defined as Mi=Y^^i''ni. Mi is proportional to the con- 
stant total mass of the system and was chosen equal 1. 
Sc and tc were determined according to Eq. ((SJ. Fig. 
clearly shows the crossover from fast power-law cluster 
growth to slow aggregation. Moreover, Eq. |(SJ) is con- 
firmed by the data collapse. In addition, we checked rate 
coefficients obeying Eq. with different values of A and 
a and found in all cases agreement with Eq. (O, as well. 

Comparison of the dimensionless rate lO with Eq. (0) 



shows that the time unit is t* 



3r, 



where h is the 



■n2feBT' 

initial number density of primary particles. As an ex- 
ample, we compute the physical time scale t* for pos- 
sible realistic parameters. With ?7=10~*Pas, T=300K, 
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FIG. 1: The average size s{t) scaled with Sc vs. time t scaled 
with tc for different values of k. The inset shows the original 
data, and the dashed line is the behavior for an uncharged 
system (fc=0). During integration the time step could be in- 
creased, since the rates for mergers become small in the stable 
regime. Thus, t/t* is larger than the number of time steps. 



a*=l/im and volume fraction i^=lQ^^ {n^2-10^mm^^) 
one obtains i*wl50s. Assuming er~l, charging with a 
single elementary charge on each primary particle corre- 
sponds to a value A:«0.2. For these parameters one finds 
with the dimensionless results of Fig. ^that the crossover 
to slow aggregation occurs within hours. 

In the unstable regime the system behaves basically 
as in the uncharged case and the cluster size distribu- 
tion is approximately given by Eq. (0). Fig. |21 displays 
a double logarithmic plot of the scaled size distribution 
in the stable regime for three different times and two dif- 
ferent choices of the rate coefficients. The coagulation 
kernel for the dashed lines obeys Eq. (O with fc=0.1 and 
a=l/3 (Brownian coagulation of charged spheres). Solid 
lines correspond to Eq. © with fc=0.05, a=0 and A=— 2. 
The inset shows the latter data unsealed (ti<i2<^3)- The 
parameter a—0 corresponds to the limiting case, where 
the Coulomb repulsion of large clusters becomes inde- 
pendent of the cluster radius. For comparison the dot- 
dashed line shows a decay ^ x~^. The labels a, b, (j){a) 
and 0(6) mark analytical expressions for the asymptotic 
edges of the distribution and the corresponding heights, 
which will be derived in the following. It can be seen that 
for monopolarly charged systems the cluster size distri- 
bution asymptotically converges to a self-similar solution 
analogous to Eq. Furthermore, Fig. El strongly sup- 
ports the claim that (j){x) is independent of the details of 
the underlying coagulation kernel. The data suggest that 
the scaling function (j){x) decays as (j){x)^x~'^ between 
sharp edges and is zero outside. As can be observed in 
the inset of Fig. [21 a consequence of this is that Ui is 
independent of the average size s between the edges of 
the size distribution. Based on the numerical data we 




i / s(t) 

FIG. 2: Double logarithmic plot of the scaled cluster size 
distribution for different times and choices of the coagulation 
kernels. The inset shows original data. 



assume 



(j)Q X for X G [a, b] 
for X ^ [a, b] 



(9) 



Applying Eq. © one finds A/o=(0o/s) (l/a-1/6) and 
Mi=l=(j>o\n{b/a), together with s—Mi/A'Iq leading to 



ln(5/a) = [1/a- 1/6] 



(10) 



The characteristic feature of the coagulation rates 
for monopolarly charged systems is that (in the stable 
regime) coagulation is exponentially suppressed with in- 
creasing cluster sizes. This feature does not depend on A, 
k, and a or whether one applies Eq. © or Eq. ©. The 
exponential suppression causes the rates for coagulation 
of clusters located at the left edge of the size distribu- 
tion to be much larger than rates for any other coagu- 
lation events. Thus, clusters at the left edge with index 
i^[as,as+Ai] coagulate, and mass conservation requires 
that this leads to an increase of the number of clusters 
with i€[2as, 2as+2Ai]. Hence, we expect that asymptot- 
ically b—2a. Between the left and the right edge of the 
size distribution n.i is essentially not changed. The scal- 
ing function Eq. (jSJ is the only possibility to guarantee 
such a behavior. Inserting b=2a into Eq. H10|l leads to 
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a= 



2 In 2 
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in good agreement with the numerical data. Fig. [3 

An important consequence of the universal form 
Eqs. (|9I11() is that the relative width of the cluster size 
distribution = ^ (MiMqIM\') — 1 approaches a univer- 
sal value (7^ with 

CTr — > = Vl/[2(ln2)2] - 1 « 0.2017 . (12) 
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FIG. 3: Relative width cr,- of the cluster size distribution for 
monopolarly charged and corresponding uncharged systems. 
k = (0.5, 0.1, 0.05) for A = (0, -2, -4). 



A=0, fc=0.2 and different values of a, in accordance with 
the predicted asymptotic behavior in Eq. (|14|l . 
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FIG. 4: s{t) times [ln(f)] ^-c for different values of a. 



Fig-Elshows numerical data for (Tr with A=(0, —2, —4) for 
monopolarly charged particles (solid lines) in comparison 
to the corresponding uncharged cases (dashed lines). (tJ?° 
is considerably smaller than the relative width for any 
physical system we found in the literature. 

For an initially narrow cluster size distribution the rel- 
ative width first grows similar to the uncharged situation 
until time tc, when it starts decreasing again, induced 
by the exponential suppression of further aggregation. 
We refer to this as self-focussing. The example discussed 
above, with t*ssl50s, shows that this self-focussing effect 
can occur within an experimentally accessible time. 

In the following, we specify what kind of time depen- 
dence the average size s exhibits in the stable regime. 
For this purpose we calculate the time derivative j^Mq = 
Si According to the previous results, in the sta- 

ble regime -^Mq reduces to a transfer from the left edge 
of the size distribution (iwas) to the right edge (iw2as). 
This behavior together with Eq. ^ leads to -^Mq = 
— ^ Rasas III tlis stable regime the coagulation 

kernels asymptotically obey Ry^y ~ exp(— fc^y^""). Us- 
ing s — Mi/Mq we arrive at 



ds{t) 
dt 
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-Ks{ty 



with K 



(13) 



This results in the time interval t2 — ti « I{s2) — ^i^i) 
in which the average cluster size increases from si to 
S2, with I{y) = dx x"^ exp{Kx^~°'). Application of 
I{y) Ri y'^ exp{Ky^~'°') for large y^l, with the exponent 
7 depending on a, leads to s(t) ~ (Ini — ■^\ns{t))^^ . 
Thus, we can conclude that asymptotically 



s{t) - (hit): 



(14) 



A similar long-time behavior was observed in [1, [ifl, IHl . 
Fig. 21 shows numerical data for Rij obeying Eq. ijHJl with 



In conclusion, we investigated Smoluchowski's coagula- 
tion equation for monopolarly charged particles. Charge 
accumulation leads to a crossover from power-law to sub- 
logarithmic cluster growth. Asymptotically the cluster 
size distribution approaches a universal scaling form, in- 
dependent of any parameter included. This implies a 
self-focussing of the relative width of the distribution. A 
first experimental check is encouraging jl2l |. 
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